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Abstract. We study the density of the roots of the derivative of the characteristic 
polynomial Z{U, z) of an N x N random unitary matrix with distribution given by 
Haar measure on the unitary group. Based on previous random matrix theory models 
of the Riemann zeta function C{s), this is expected to be an accurate description for 
the horizontal distribution of the zeros of C'(s) to the right of the critical line. We 
show that as ^ oo the fraction of the roots of Z'{U,z) that lie in the region 
1 — x/{N — 1) < |z| < 1 tends to a limit function. We derive asymptotic expressions 
for this function in the limits x —> oo and x — > and compare them with numerical 
experiments. 
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1. Introduction 

We study the density of the roots of 

d d ^ 

Z'(U, z) = — det( Iz ~U) = — Y\ (z- e'^0 , z e C, 
az dz 

i=i 

where f/ is a random NxN unitary matrix, with respect to the circular unitary ensemble 
(CUE) of random matrix theory (RMT). Our main motivation is to investigate the 
horizontal distribution of the zeros of the derivative of the Riemann zeta function. 
The zeta function is defined by 

oo ^ 

^ — ^ n'^ 

n=l 

and has an analytic continuation in the rest of the complex plane except for a simple 
pole at s = 1. There are infinitely many non-trivial solutions to the equation ({s) = 
in the strip < a < 1; the Riemann hypothesis (RH) states that they all lie on the 
critical line a = 1/2. The interest in the horizontal distribution of the zeros of C'{s) is 
motivated by its connection with RH. In 1934 Speiser [|I| showed that RH is equivalent 
to C'{s) having no zeros in the region < a < 1/2. Furthermore, up to now the most 
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efficient ways of computing the fraction the of zeros of the Riemann zeta function on the 
the critical hne are based on what is known as Levinson's method 0; it turns out that 
the zeros of C'{s) close to the critical line have a significant effect on the efficiency of 
this technique , therefore it is important to know how they are distributed. Levinson 
and Montgomery ^ proved a quantitative refinement of Speiser's theorem, namely that 
and C'{s) have essentially the same number of zeros to the left of a = 1/2, and 
showed that as T ^ oo, where T is the height on the critical line, a positive proportion 
of the zeros of C'{s) are in the region 

^ 1 , NloglogT 
2 logT 

Subsequent improvements of Levinson and Montgomery's results, first by Conrey and 
Ghosh PI, then by Guo p, Soundararajan and recently by Zhang [0] have established 
that a typical zero of C'{s) tends to be much closer to the critical line and that 
conditionally on RH a positive proportion lie in the region 
1 C 



2 logT' 

for some positive constant C . Their distribution, however, is still unknown. Other 
results on the zeros of C'i.^) can be found in [||. 

Over the past thirty years, overwhelming evidence has been accumulated which 
suggests that the local correlations of the non-trivial zeros of coincide, as T ^ oo, 
with those of the eigenvalues of hermitian matrices of large dimensions from the Gaussian 
unitary ensemble (GUE) |Q. As ^ oo, the GUE statistics are in turn the same as 
those of the phases of the eigenvalues of N x N unitary matrices, on the scale of their 
mean distance 27r/iV, averaged over the CUE ensemble. More recently, however, it was 
realized that RMT not only describes with high accuracy the distribution of the Riemann 
zeros, but that it also provides techniques to make predictions and computations about 
the Riemann zeta function and certain classes of i^functions that previous methods 
had not been able to tackle. This started with the work of Keating and Snaith on 
moments of the Riemann zeta function and other L-functions. Their key observation 
was that the locally determined statistical properties of ({s) high up the critical line 
can be modelled by characteristic polynomials Z{U, z) of random unitary matrices U . 
In this model the two asymptotic parameters, T for C,{s) and for f/, are compared by 
setting the densities of the zeros of and of the eigenvalues of U equal, i.e. 

T 

N = log — . 



This approach has since been extremely successful [11 



Following the same ideas, in this paper we suggest that the density p{z) of the roots 
of Z'{U,z) will accurately describe the distribution of the zeros of ('{s). A classical 
theorem in complex analysis states that if p{z) is a polynomial, then the roots of p'{z) 
that are not roots of p{z) lie all in the interior or on the boundary of the smallest convex 



polygon containing the zeros of p{z) (see, e.g., [T^). Therefore, since the eigenvalues of 
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a unitary matrix have modulus one, the solutions of the equation Z'{U, z) = that are 
not zeros of Z{U, z) are all inside the unit circle. If s = 1/2 + it + u, t G M, denotes the 
point at which ('{s) is evaluated, then the region of C to the right of the critical line 
is mapped inside the unit circle by the conformal mapping z = e~". Thus, the radial 
density 

\z\ p{z) d(f) 





becomes the analogue of the horizontal distribution of the zeros of C'{s) to the right 
of the line a = 1/2. However, instead of p{z), it turns out to be more convenient to 
consider Ip(x), the fraction of the roots in the annulus 1 — x/{N — 1) < \z\ < 1, where 
X is the scaled distance from the unit circle. Our main results concern the asymptotics 
of Ip(a;). We show that as increases the roots of Z'{U,z) approach the unit circle, 
and Ip(x) tends to a function independent of N. Furthermore, we obtain the following 
asymptotics as — oo: 

Ip(x) ~ 1 — 1/x, X — > oo, X = o{N), 

^ 97r 22577 22057r ^ ^ 

These formulae are then tested numerically. 

The paper is organized as follows: in section Q we introduce the mathematical 
problem and describe the main properties of the roots of Z'{U,z); in section ^ the 
asymptotics of Ip(x) as a; ^ 00 are computed; using heurstic arguments, in section |, 
lp(a;) is derived in the limit x — > 0; section |^ concludes the paper with final remarks. 

2. The distribution of the roots of Z'{U,z) 

The CUE ensemble of RMT is defined as the space U(A^) of N x N unitary matrices 
endowed with a probability measure dfi{U) invariant under any inner automorphism 

U ^ VUW, U,V,W e V{N) . 

In other words, d/i(?7) must be invariant under left and right multiplication by elements 
of U(A^), so that each matrix in the ensemble is equally weighted. There exists a unique 
measure on the unitary group U(A^) with this property, known as Haar measure. The 
infinitesimal volume element of the CUE ensemble occupied by those matrices whose 
eigenvalues have phases lying between = {9i,92, ■ ■ ■ , 9n) and 6 + d'^d is given by 

n{N)A\o)d^e, (2.1) 



where 



A{e)= n and m) = 7A 



<k<N 

Our goal in this paper is to study the density of the roots of Z'{U,z), where U is a 
random unitary matrix with distribution given by (p.l|) . The analogous problem for the 
Ginibre ensemble has been studied by Dennis and Hannay |jl3 . 
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(a) iV = 20 (h) N = 50 

Figure 1. Zeros of characteristic polynomials of random unitary matrices (□) and of 
their derivatives (a). 



In figure |l] are plotted the zeros of tlie cliaracteristic polynomials Z{U, z) and 
of their derivatives of two unitary matrices taken at random with respect to Haar 
measure for = 20 and = 50. Such matrices can be easily generated numerically 
by taking complex matrices whose entries are independent complex random numbers 
with Gaussian distribution, and then by applying Gram-Schmidt orthogonalization to 
the rows or columns (see, e.g., fT^). There are a few qualitative features that can be 
immediately observed. Firstly, since the distribution ( ^.\\) is translation invariant on the 
unit circle, the density of the roots of Z'{U,z) depends only on the distance from the 
origin. Secondly, as mentioned in the introduction, the roots of Z'{U,z) are all inside 
the unit circle. This property can be easily understood with the following argument. 
Let zi, Z2, ■ ■ ■ , zn be complex numbers; if they all are on the same side of a straight 
line passing through the origin, then 

11 1 

zi + Z2 + ■ ■ ■ + zn 7^ and — + — + ••• + — ^0. (2.2) 

Zi Z2 Zn 

Let {zj}^^i be the roots of a polynomial p{z) and z a point outside the smallest convex 
polygon containing {zj}^^^. Now, consider a straight line passing through z and lying 
outside such a polygon. Because of equation (|2.2|), the logarithmic derivative 



p'iz) 1 1 1 



p{z) Z — Zi Z — Z2 Z — Zn 

cannot vanish. There are two others less obvious features that figure |I| reveals and that 
become more apparent as N increases: firstly, the roots of Z'{U,z) concentrate in a 
small region in proximity of the the unit circle; secondly, given two consecutive zeros 
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of Z{U,z), say e^^^ and e^^^+i, close to each other, there often appears to be a root of 
Z '{U, z) near the midpoint 



2 - >2.3) 

In the following sections we shall give quantitative interpretations of these properties. 

Let us set u = Rez and v = lm.z] furthermore, denote by {\j{6)}^~l the set of 
roots of Z '{U, z) and consider the linear functional 

F[r],^^e) = [ r{z)6 {z - X,{0)) d'z = r {X,{0)) , 

where (Pz = du dv and t{z) is an infinitely different iable complex function whose partial 
derivatives with respect to u and v decrease faster than any power of 1/ \z\. Moreover, 

5{z- X,{0)) =6{u-Re (A,(0))) 6{v-lm (A,(0))) 

is the product of two Dirac delta functions with real arguments. Then, we have 



T[z)p{z)d'^z 



c 



iV- 

where the distribution 



S^E / / r{z)6{z-X,{e))A\e)d^zd^9, 



p(') ■■=jr^Il [ ^^■(^)) ^'(^)d''^ (2-4) 

defines the density of {Xj{6)}^S'i . 

The main tool that we shall use to evaluate ( p.4| ) is a basic identity that expresses 
Toeplitz determinants in terms of integrals over the unitary group. If 

oo 
fc=— oo 

is a complex function on the unit circle, then we denote by Dn-iIJ] the determinant of 

the Toeplitz matrix 

/ /o /i • • • fw-i \ 

\ f-{N~l) f-{N-2) ■ ■ ■ fo / 

Now, let $j be a class function, i.e. a complex function on U(A^) such that 

$/ {VUV-^) = $/(f/), V,U e V{N). 
Furthermore, suppose that 

$,Kf/) = /(^i)/(^2)---/(M, (2.5) 
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where {e'^-' }jLi are the eigenvalues of U. The Heine-Szego identity |T5[ states that 



Q{N) 



[2.6) 



In order to apply this formula, we must express the sum of delta functions in ( |2.4| ) 
as a product of the form ( |2.5| ). The zeros of Z'{U,z) that are not multiple roots of 
Z{U, z) are the same as those of the logarithmic derivative of Z{U, z). Since the set of 
unitary matrices with degenerate eigenvalues has zero measure, we rewrite the integrand 
in equation (|2.4|) as 

N-l 



Y,6{z-X,) = d{Z'{U,z)/Z{U,z)) 



dz 



[Z'{U,z)/Z{U,z)] 



[2.7) 



In the next step we use the integral representation of a delta function: 



The complex delta function in the right-hand side of equation ( p.7|) now becomes 

'i (z'iU,z] 



6{Z'{U,z)/Z{U,z)) 



47r2 



exp 



ZiU,z) 



w 



Z'{U,z) 
Z{U,z) 



w 



di'w. 



(2. 



Clearly, the identity (|2.6|) can be applied to the argument of the integral (|2.8|) ; 
furthermore, the Jacobian in equation (|2.7| ) can be transformed into a product of the 
form (|2.5| ) by using the following representation of the modulus square of a complex 
number: 



da' 



G{a, z) 



a G 



a=0 



where 



G{z, a) := exp [ia {z + z) /2] + exp [a {z — z) /2] 
Finally, the density (|2.4| ) becomes 

Piz) - 



X 



47r2 {N - 1) 

J (^DN-i[exp{ig)]{w,a,z) + DN-i[exp{ih)]{w,a,z)^d'^w 



[2.9) 



a=0 



where 



and 



g{e;w,a,z) := - 



1 



w 



w 



z -e'^ z - e^i" 



h{9; w, a, z) := - 
are real functions. 



w 



+ 



w 



2 \z — e'^ z — e-i^ 



a 
2 



+ 



1 



z-e'^)^ (z-e-'^)' 
1 1 
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Computing ( |2.9| ) exactly appears to be a formidable task. However, there exists a 
powerful theorem of Szego that will allow us to compute the leading-order asymptotics 
of p{z) as — * oo. 

The strong Szego limit theorem. Let 



k= — oo 



be a complex function on the unit circle. If the series 



k=—oo 



k=—oo 



converge, then 



DN^i[exp{r])] = exp fjoN + fj.kfjk + o(l) 



N ^ oo. 



(3.1) 



(3.2) 



k=l 



The first proof of this theorem was given by Szego in |TB[ under stronger conditions; 
several proofs have since been developed fl^, [Tsll - 

In the region \z\ < 1, g{9; w, a, z) and h{6; w, a, z) are just sums of geometric series 
and of their derivatives whose Fourier coefficients can be easily computed. We have 



g{9;w,a,z) 



J2 h^"^^ and h{e- w,a, z) 



k=—oo 



k = — OD 



fee 



ike 



where 



and 



9o = 0, gk = -—z 



^^k-l iA^fe-2 



A; G 



ho = 0, hk = -—z 



^ (A; -1)^^-2, A;gZ+. 



Since g{9;w,a,z) and h{9;w,a,z) are real, g^k = dk and h_k = ^fe- Computing the 
argument in the exponential of equation ( p.2|) involves only summing and differentiating 
geometric series. We obtain 

a Re(ryz) 



9k\ 



E{g) = J2k 

k=l 

oo 

E{h) = J2k\h 



k=l 



where r = |z| and 

c{w, a, r) 



\w\ 



4 (l-r2) 



c{w, a, r) + 



c{w, a, r) + 



a 



a Im(wz) 



(3.3a) 
(3.35) 



3r^ 



+ 



2^3 
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10 12 14 16 18 20 



10 12 14 16 18 20 



(a) 



(b) 



Figure 2. Comparison (a) and difference (b) between Ip(a;) computed numerically for 
N = 800 (0) and formula (U) ( ). 



As a consequence, the second sum in ( |3Tl| ) is finite. Furthermore, we have 



hm 


gk+1 


= hm 


hk+1 


= r < I 


fe— >oo 


9k 


k^oo 


hk 





Hence, the first series in ( |3.1| ) converges too, and the strong Szego hmit theorem apphes. 

When computing derivatives of the asymptotics of the integral (|2^ ), care must be 
taken that the error term does not become comparable or even greater than the leading- 
order term. This could happen, for example, if the remainder were a highly oscillatory 
function of a. It turns out that the convergence of DN-i[exp{ig)] and Dj^_i[exp{ih)] 
to ( p.2|) is so fast that the derivatives of the error term remain small. This is proved in 
the appendix. 

In equation ( p.9|) the second derivative commutes with the integral, hence 
differentiating twice with respect to a gives 

r 

^-^ Dn-1 [exp(i5f)] (w, a, z) + Dn-i [exp(i/i)] (ly, a, z) 



I |2 2 

\w\ r 



6r^ 



{l-r^f {l-r^f {l~r^f 



a=0 



exp 



\w\ 



4fl_^2^2 



This expression can be trivially integrated. Finally, we obtain 

2 



\2 ' 



N 



oo. 



7r(A^- 1) {l-r^y 

As anticipated, p{z) depends only on the distance from the origin r and is asymptotically 
concentrated in a small region near the unit circle, which explains the migration of the 



Random matrix theory and the zeros of ('{s) 



9 



roots of Z '{U, A) observed in figure as increases. Since p{z) is a density, it must be 
normalized to one, tfierefore we require 

rp{z) drd(f) ~ 1, oo, e 0. 



Let us now define Ip(a;) to be the fraction of tfie zeros in tlie annulus 
1 — x/ [N — 1) < r < 1, where x = o{N), i.e. 



Ip(x) 



2tt />1-x/(N-1) 



rp{z) drd0 ~ 1 
1/x, N ^ oo, X oo. 



N 



-1 l-x/{N-l) 



As ^ oo the leading-order term of Ip(x) is independent of A^. Equation 



(3.4) 
is the 



main result of this section. In figure Q formula 
for a matrix of dimension A^ = 800. 



is compared with Ip(x) computed 



4. The asymptotics x ^ 

The small x asymptotics of Ip(a;) requires first the evaluation of the limit x — and then 
of the limit A^ oo. Szego's theorem gives an asymptotic expression as A^ — * oo before 
the limit x —>■ can be taken, therefore provides information only for relatively large x. 
In order to determine Ip(x) in the limit x — 0, the integral in equation (|2.9|) needs to be 
evaluated for finite A^. Such computation seems to be an extremely difficult task: the 
integrand has essential singularities that usual techniques in RMT and complex analysis 
cannot tackle. Notwithstanding such obstacles, it turns out that Ip(x) can be derived 
in the limit x — with the help of a heuristic argument. 

As was mentioned in section |^, from figure |1| it appears that if e^^^ and e'^-'+^ are 
two consecutive roots of Z{U, z) which are close to each other, then as A^ — * oo there 
is often a root of Z'{U,z) near the midpoint ( p.3| ). Thus, one might be led to assume 
that for small x the distance from the origin r is distributed like 



x/(A^ 



gie, + e^e.+i ^ ^2 + 2 cos {27iS/N) 



2 



S = —\dj+i-dj\, 1<J<N. 



o , ^< J < N, (4.1) 

where S is the rescaled distance, or spacing, between phases of consecutive eigenvalues, 
i.e. 

A^ 
2^ 

This is trivially true only for N = 2; since S = 0(1), for N > 2 equation ( [4.1|) would 
imply an average distance of the zeros of Z'{U,z) from the unit circle of order l/N"^, 
and therefore a dependence of Ip(x) on A^ even at the leading order, which contradicts 
the numerics reported in figure ^ and formula ( p.4| ). 

It turns out that behaviour of Ip(x) as x — can be understood using Dyson's 
electrostatic model for the CUE ensemble (see, e.g., [|l9l). The zeros of Z{U,z) can be 
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Figure 3. Schematic representation of the two components £^((7) and Eab(<?) of the 
electric field (O) at a point q close to the unit circle in Dyson's electrostatic model. 



interpreted as N unit charges confined in a two-dimensional universe and moving in a 
thin circular conducting wire of radius one. The electric field generated at any point in 
the complex plane by this Coulomb gas is just the complex conjugate of the logarithmic 
derivative of Z{U,z), i.e. 

E(z) = ^ + ^ + 



z — e"^^i z — e~'^2 z — e-i^iv ' 

Hence, the zeros of Z'{U,z) are located where E(^) vanishes. Now, the field at point 
q, whose distance from the unit circle is of order can be separated into two 

components: the first one is the field of the two charges closest to q, let us denote 
them by A and B, whose strength is clearly of order N; the second one is the field 
generated by the other N — 2 charges. As N increases, q approaches the unit circle, and 
the latter component of E(g) can be approximated by the field of a continuous circular 
charge distribution with density /i = N/{27i), i.e. 

E(g)^E^(g) + EAB(g), (4.2) 

where E^(g) is the field generated by fi and Eab(9) the one determined by A and B. For 
large N, E^(g) and Eab(?) have approximately opposite directions, thus E(z) might 
vanish at q only if |E^(q')| = O(A^). This situation is described schematically in figure 0. 
Determining the order of magnitude of E^(g) is a simple exercise in electrostatics. If the 
continuous charge distribution filled the whole unit circle, the field inside it would be 
zero. Thus, by linearity E^(g) is equal and opposite to the field of a circular arc [—9, 9] 
with charge density /i and containing no more than four eigenvalues. We have 



/ M r 1 , sin^^ 



Random matrix theory and the zeros of ('{s) 



11 



where d{9) is the distance between q and e^^. By setting d{9) = t{9)/N, with t{9) = 0(1) 
and applying the mean value theorem we obtain 



Since the length of the interval [—9, 9] is of the order of few level spacings, 9 = 0(1/ N) 
and |E^(g)| =0(iV). 

The field Eab(-2) vanishes at the midpoint (|2.3D, but for large the presence of 
E^(z) shifts the zeros of E(z) at a distance of order 1/A^ from the unit circle. Since 
at first approximation the contribution of E^(z) does not depend on the spacing S, 
from (|4.1|) it is reasonable to assume that as a; — ^ the distance of the zeros of Z '{U, z) 
from the unit circle is approximately distributed like S'^/N. Hence, we shall conjecture 
that 

^-13^. a;^0, (4.3) 

where /3 is a constant independent of A^. 

From equation (|4.3|) it is straightforward to derive an asymptotics expression for 
Ip(a;) as X — > 0. We have 



Ip(x) ~ / 7(?/) dy, a; ^ 0, 

^0 



where 



1/2- 



and -Pcue(5') is the CUE spacing distribution in the limit N ^ oo, i.e. 

PcMS) := lim PcuE(iV;^). 
-Pcue('S') has a power series expansion with infinite radius of convergence: 

oo 

PcMS) = ^(/ + 4)(/ + 3)Ei S'+'. (4.4) 

(=0 

There exist efficient algorithms for computing the coefficients Ei (see, e.g., |2^), which 
with symbolic mathematical packages can be evaluated exactly up to very high values 
of I. Using (|4.4|) one can easily obtain a series expansion for Ip(x): 

oo / 2 \^ 

IP(^)~5Z(^) {l + ^)Eix'-^, x^O. (4.5) 

Now we have to determine the parameter f3, which can be found only empirically. It 
turns out that if we set /3 = 1/2 there exists an astonishing agreement between ( |4.5| ) 
and numerics in the region where the large x asymptotics is not valid. This is shown in 
figure |. 

Notwithstanding the accuracy with which the series ( [4.5|) models the numerical 
data, it can be expected to approximate Ip(a;) only for small x; indeed, it tends to one 
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much faster than 1 — 1/x. Furthermore, in deriving ([4 .51 ) we have imphcitly assumed 
that the main contribution to Ip(x) comes only from the two zeros of Z{U,z) closest 
to a given root of Z'{U, z); we now need to estimate in what region this assumption is 
justified. 

Let us consider k + 2 successive eigenvalues of a unitary matrix. Since Haar measure 
is translation invariant on the unit circle, the corrections to (|4.5| ) will depend by all 
possible rescaled distances 
N , 

Sk = — \0j+k+i-0j\, l<j<N. 

271 

However, as 5*^ — > 0, the limit densities Pc\jE{Sk) go to zero very fast, hence the 
dependence on Sk does not affect the first few terms of the series ( |4.5| ). For example, 
consider three consecutive zeros of Z{U, z) close to each other, say 

exp , exp (a + S^)^ and exp {a + Si)^ , Si>Sq>Q. 

Simple algebra and Dyson's model indicate that the rescaled distance x should be 
distributed like 

/3'(5o + 

for some constant [3' independent of A^. It turns out that 

PcuE(^i) = ^ + 0(5f), (4.6) 



which suggests that the contributions to ( |4.5|) due to 5*1 leave the coefficients of the 
terms x'^ with d < A unchanged. For > 1, PcvE^Sk) goes to zero as 5*^ ^ even faster 
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(a) 



(b) 



Figure 5. (a) Fraction of zeros of Z'{U, z) in tlie region 1 — x/{N — 1) < |z| < 1 for 
N = 800 (0), N = 500 (O ) and = 400 (a). (b) Same data as in (a) for N = 800 

(0) compared with formula (3.4) ( ) and with the series (4.5) truncated at Z = 30 

and with /3 = 1/2 ( ). 



than ( [4.6|) . These considerations and formula ( [4.5|) give the following expression for the 
asymptotic expansion of Ip(x): 

64 128 



Ip(x) = — x^/^ 



X 



X 



7/2 



+ 0(x^). 



2257r 22057r 

However, from the agreement with numerics observed in figure ^ and p(b)| , we would 
expect that the corrections to the series ( |4.5| ) should be negligible up to terms of order 
much higher than x'^. 

The integrated distribution Ip(x) is plotted in figure |5(a)| for zeros of Z'ilJ^z) 
computed numerically for random unitary matrices of various dimensions; clearly, Ip(x) 
tends to a limit function. Figure |5(b)| shows that equation ( |3.4| ) and the series ( |4.5| ) 
together, although asymptotic formulae, approximate Ip(x) with high accuracy for all 
X > 0. 



5. Concluding remarks 

We have investigated the density p{z) of the roots of Z'{\J,z), where Z{\J,z) is the 
characteristic polynomial of a random N x N unitary matrix with distribution given by 
Haar measure on the unitary group. Since the locally determined statistical properties 
of the Riemann zeta function high up the critical line can be modelled by Z{U,z), it 
is expected that p{z) will accurately describe the behaviour of the zeros of C'i^)- It 
turns out that instead of p{z), it is more convenient to study Ip(a;), the fraction of the 
roots in the region 1 — x/{N — 1) < |z| < 1. In the analogous problem for the zeta 
function, this is equivalent to looking at the fraction of the zeros of C'{s) in the region 
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1/2 < a < 1/2 + x/ logT where a — Re(s) and T is the height on the critical hne. It is 
shown that as N oo, Ip(,x) becomes independent of N. 

The density p{z) can be defined as the average over U(A'") of the sum of delta 
functions 

where the Xj{0) are the zeros of Z'{U,z). The behaviour of Ip(a;) for large x can be 
computed by applying standard techniques for integrals over U(A'^). The sum of Dirac 
deltas can be manipulated in such a way that eventually the average over U(A^) is 
reduced to the computation of the second derivative of an integral over the complex 
plane of the sum of two Toeplitz determinants. Furthermore, the integrand satisfies the 
hypothesis of the strong Szego limit theorem, which gives a simple asympotic expression 
for such determinants. Further simple manipulations lead to 

lp(a;) ~ 1 — l/x, X = o{N), x oo. 

The limits a; ^ and ^ oo do not commute; the evaluation of the asymptotics 
of lp(a;) as X ^ requires first the computation of the limit a; — ^ and then of the limit 
oo. The application of Szego's theorem clearly prevents this, therefore provides 
information only for relative large x. However, Dyson's electrostatic model for the CUE 
ensemble leads naturally to the assumptions that for small x, as N ^ oo, the distance 
of the roots oi Z'{U,z) from the unit circle is on average of order 1/N and is distributed 
like the square of the spacings between phases of consecutive eigenvalues of unitary 
matrices in the CUE ensemble (appropriately rescaled). These two simple hypotheses 
give a conjecture for Ip(a;) as x — > whose agreement with numerical experiments covers 
with high accuracy the region where the large x asymptotics fails. 

Unfortunately, the zeros of ('{s) are poorly understood, and there is not even a 
conjecture for their horizontal distribution to compare with the results derived in this 
paper. Given that it seems extremely difficult to obtain an analytical expression of such 
a quantity, we believe that it would be interesting and worthwhile to conduct a thorough 
numerical study as independent verification of the model presented here. 
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In this appendix we show that in the cases considered in the present paper, the error 
term in formula ( |3.2| ) remains small when differentiated with respect to a; in other 
words, we want to prove that 



DN_i[exp{ig)]{a) + DN-i[exp{ih)]{a) 



a=0 



exp(-E((7))+exp(-^(/i)) 



oo, 



(Al) 



J a=0 



where E{g) and E{h) are defined in equations (3.3)||. The main idea of the proof is 
quite simple: we first represent the Toeplitz determinants with exact formulae and 
differentiate them with respect to a; then we take the limit N ^ oo. We shall consider 
only DAr_i[exp(i(7)], as the proof for D ^ _i[exp{ih)] is completely analogous. 
The identity (|2.6| ) allows us to write 



D 



N-l 



[exp{ig)]{a) = [ exp (i V ^^(a) Tr(t/'=) ) dfi{U). 

-^UW V k=-oo J 



Replacing the exponential function by its power series gives 



CO CO 



(i^fc(«)Tr([/'=)) 



bk=0 



-igk{a) Tr([/'=))'"= 



dfi(U). 



bk=0 

Integrals over U(A) of product of traces of unitary matrices have been computed by 
Diaconis and Shahshahani [^. Let Xj be nonnegative integers such that 

Al > A2 > . . . > A, 

and 

L = Al + A2 H h As = lai + 2a2 H h ra^, 

where denotes the number of times the integer k appears among the Xj s. We call 
Xa = (Al, A2, . . . , As) a partition of L. Consider the integral 



n (Tr(f/*^))"'= (Tr(f/*^))'^- d/i([/). 

U(^) k=l 



HXaAb) 

It turns out that /(A^, Xb) = unless A^ = A;,; furthermore, we have (see, e.g., 

r 

I{Xa,Xb) = Sx,^x,l[k^''ak\ if L<N, 



k=l 
r 



(A2) 



/(A., A,) < 5,^,, n 



if L > A. 



A;=l 



I To simplify the notation and emphasize the dependence on a, we omit the variables w and z in this 
appendix. 
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Therefore, we obtain 

D^^,[exp{ig)]{a) = J^H^a, Xa)U(-ir ( i ^ (^3) 

where the sum is over all partitions and the product over all the integers (without 
multiplicity) of a given partition. Because of equation (^), asymptotically formula 
tends to 

\a k A;=lafc=0 

OD 

2 



exp - ^ A; \gk 



a) 



k=l 



This proof of the strong Szego limit theorem was derived by Bump and Diaconis [|1^ . 

In order to prove (^) we need to take the second derivative with respect to a 
of (|A3| ) and then the limit N oo. Differentiating the right-hand side of equation ( [^ ) 
is tedious but elementary. We shall carry out only the first derivative, since the second 
one is completely analogous. We have 

^-2^/(A.,A.)^(-l) 5, (a.-l)!a.! 



da 



j^k '^^r) 
In the limit N ^ oo the right-hand side of ( [A4|) becomes 

da I ~Z^^ )= d^^^^^'^^''' 

fc=l V k=l J 

which is the same expression obtained by differentiating Szego's strong asymptotic 
formula. Similarly, differentiating (^5) and then taking the limit N oo gives 



d'^DN-i[exp{ig)]{a) 



da 



2 



fdE{g)y d'E{g) 



\ da J da^ 
X exp {-E{g)) , N ^ oo. (A5) 

This completes the proof of (|A1|) . 
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